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Binary collisions between charged particles in an external magnetic field are considered in second- 
order perturbation theory, starting from the unperturbed helical motion of the particles. The 
calculations are done with the help of an improved binary collisions treatment which is valid for any 
strength of the magnetic field, where the second-order energy and velocity transfers are represented 
in Fourier space for arbitrary interaction potentials. The energy transfer is explicitly calculated for 
a regularized and screened potential which is both of finite range and non-singular at the origin, and 
which involves as limiting cases the Debye (i.e., screened) and Coulomb potential. Two distinct cases 
are considered in detail, (i) The collision of two identical (e.g., electron-electron) particles; (ii) and 
the collision between a magnetized electron and an uniformly moving heavy ion. The energy transfer 
involves all harmonics of the electron cyclotron motion. The validity of the perturbation treatment 
is evaluated by comparing with classical trajectory Monte-Carlo calculations which also allows to 
investigate the strong collisions with large energy and velocity transfer at low velocities. For large 
initial velocities on the other hand, only small velocity transfers occur. There the non-perturbative 
numerical classical trajectory Monte-Carlo results agree excellently with the predictions of the 
perturbative treatment. 

PACS numbers: 03.65.Nk, 34.50.Bw, 52.20.Hv, 52.40.Mj 



I. INTRODUCTION 

In the presence of an external magnetic field B the problem of two charged particles cannot be solved in a closed 
form as the relative motion and the motion of the center of mass are coupled to each other. Therefore no theory 
exists for a solution of this problem that is uniformly valid for any strength of the magnetic field and the Coulomb 
force between the particles. The energy loss of ion beams and the related processes in a magnetized plasmas which 
are important in many areas of physics such as transport, heating, magnetic confinement of thermonuclear plasmas 
and astrophysics are examples ofphysical situations where this problem arises. Recent applications are the cooling 
of heavy ion beams by electrons [J |3, [1, S @ and the energy transfer for heavy-ion inertial confinement fusion (ICF) 
(see, e.g., for an overview). The classical limit of a hydrogen or Rydberg atom in a strong magnetic field also falls 
in this category (see, e.g., 0] and references therein) but in contrast to the free- free transitions (scattering) the total 
energy is negative here. 

Numerical calculations have been performed for binary collisions (BC) between magnetized electrons [1, and for 
collisions between magnetized electrons and ions [13, [H m, [H, In general the total energy E of the particles 
interacting in a magnetic field is conserved but the relative and center of mass energies are not conserved separately. 
In addition, the presence of the magnetic field breaks the rotational symmetry of the system and as a consequence 
only the component of the angular momentum L parallel to the magnetic field L|| is a constant of motion. So, the 
constants of motion E and L\\ reduce the phase space of the relative motion. A different situation arises for the BC 
between an electron and uniformly moving heavy ion. As an ion is much heavier than an electron, its uniform motion 
is only weakly perturbed by collisions with the electrons and the magnetic field. In this case L\\ is not conserved but 
there exists a conserved generalized energy K (T^ [T3 | involving the energy of relative motion and a magnetic term. 
The apparently simple problem of charged particle interaction in a magnetic field is in fact a problem of considerable 
complexity and the additional degree of freedom of the cyclotron orbital motion produces a chaotic system with two 
degrees (or one degree for heavy ions) of freedom [3, [H, [H, [I3| ■ 

In this paper we consider the BC between two charged particles treating the interaction (Coulomb) as a perturbation 
to their helical motions. For electron-heavy ion collisions this has been done previo usly in first order in the ion charge 
Z and for an ion at rest [l^ and in up to 0{Z^) for uniformly moving heavy ion [Illll. In Ref. [H three regimes 
are identified, depending on the relative size of the parameters a (the cyclotron radius), s (the distance of the closest 
approach), and d (the pitch of the helix). In earlier kinetic approaches [J, S S i, i only two regimes have been 



'Permanent address: Institute of Radiophysics and Electronics, 378410 Ashtarak, Armenia; Electronic address: |hrachya@irphe.am| 
t Electronic address: guenter.zwicknagel@physik.uni-erlangen.de 



2 



distinguished: Fast collisions for s < a, where the Coulomb interaction is dominant and adiabatic collisions for s > a, 
where the magnetic field is important, as the electron performs many gyrations during the collision with the ion. The 
change AEi of the energy of the ion has been related to the square of the momentum transfer Ap, which has been 
calculated up to 0{Z). This is somewhat unsatisfactory, as there is another O(Z^) contribution to AEi, in which the 
second-order momentum transfer enters linearly. Moreover, for applications in plasma physics (e.g., for calculation 
of the ion energy loss in a magnetized plasma) one calculates the angular averaged energy transfer which vanishes 
within first-order perturbation theory due to symmetry reasons and the ion energy change receives contribution only 
from higher orders ^20] . Indeed, the transport phenomena, etc., are of order O(Z^) in the ion charge. 

Here we focus on BC between two magnetized identical particles (e.g. electrons) within the second order per- 
turbation theory and its comparison with classical trajectory Monte-Carlo (CTMC) simulations. The present work 
considerably extends our earlier studies in Refs. [3, Ho] where the second-order energy transfer for an ion-electron 
collision was calculated with the help of an improved BC treatment which is valid for any strength of the magnetic 
field and does not require the specification of the interaction potential. In addition, we consider here the impact 
parameter-averaged energy transfer for the BC between magnetized electron and heavy ion moving uniformly along 
the magnetic field which has not been yet considered in Refs. [3, H^], and will give new analytical expressions which 
are more appropriate for an explicit calculation of the energy loss. Physically these two distinct cases are similar 
except the time-dependent center of mass cyclotron motion in the case of two identical particles. The paper is orga- 
nized as follows. In Sec. |TT] starting from the exact equation of motion of two charged particles moving in a magnetic 
field we discuss some basic results of the exact BC treatment for the energy and velocity transfers as well as the 
energy conservation. In the following Sec. IIIIl we discuss the velocity and energy transfer during BC of magnetized 
particles for arbitrary magnetic fields and strengths of the two-particle interaction potential. The equations of motion 
are solved in a perturbative manner up to the second order in interaction force starting from the unperturbed helical 
motion of the particles in a magnetic field. Then in Sec. IIVI we turn to the explicit calculation of the second order 
energy transfer for electron-electron collision but without any restriction on the magnetic field B. The obtained energy 
transfer involves all cyclotron harmonics of the electron helical motion. For further applications (e.g., in cooling of ion 
beams, transport phenomena in magnetized plasmas) we consider the regularized and screened interaction potential 
which is both of finite range and less singular than the Coulomb interaction at the origin and as the limiting cases 
involves the Debye (i.e., screened) and Coulomb potentials. An exact solution for two particle collision in the presence 
of infinitely strong magnetic field is considered in Sec.|Vl This also suggests an improved perturbative treatment for 
repulsive interaction and in the case of strong magnetic field. In Sec. lVII the results of the perturbative binary collision 
model are compared with CTMC simulations in which the scattering of ensembles of magnetized electrons are treated 
exactly. This comparison allows to determine clearly the range of validity of the perturbative treatment. The results 
are summarized and discussed in Sec. IVIII Some details of the calculation of the generalized Coulomb logarithm is 
described in the Appendix |^ The small velocity limits of the energy transfers are derived in the Appendix [BJ 

II. BINARY COLLISION FORMULATION. GENERAL TREATMENT 

For our description of binary collision (BC) we start with considering the equations of motion for two charged 
particles moving in a homogeneous magnetic field and the related conservation laws, in general. This will then be 
specified to the two particular case on with we will focus in this paper. The BC between two electrons and between 
an electron and an uniformly moving heavy ion. Next the quantities of interest, the velocity transfer and the energy 
transfer of particles during the binary collision, will be introduced and discussed, before we turn to the solution of 
the equations of motion in the subsequent section. 

A. Relative and cm motion and conservation laws 

We consider two point charges with masses toi, m2 and charges gie, q2e, respectively, moving in a homogeneous 
magnetic field B = Bh and interacting with the potential qiq2^^U{r) with = e^/47reo. Here Eq is the permittivity 
of the vacuum and r = ri — r2 is the relative coordinate of the colliding particles. For two isolated charged particles 
this interaction is given by the Coulomb potential, i.e. Uc(j) — 1/r. In plasma applications the infinite range of this 
potential is modified by the screening. Then U may be modeled by Ux>{v) = e~^^'^/r with a screening length A which 
can be chosen as the Debye screening length Ad, see, for example 21]. The quantum uncertainty principle prevents 
particles (for qiq2 < 0) from falling into the center of these potentials. In a classical pict ure this can be achieved by 
regularization of U{r) at the origin, taking for example C/R(r) = (l — e^*"/^) e~'^l^ jr, |2^.[23j. Here A is a parameter, 
which is usually related to the de Broglie wavelength. 
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In the presence of an external magnetic field, the Lagrangian and the corresponding equations of motion cannot 
be separated into parts describing the relative motion [r = ri — r2,v = r] and the motion of the center of mass 
(cm) [R = (miTi +m2V2)/(rni +m2),V = R], in general (see, e.g., [1, [13, [3 US] ) • Introducing the reduced mass 
= I/toi + l/m2 the equations of motion are 

v(0 - A. + 5) [v(^) X b] ^ eB - ^) [V(0 x b] + (r(t)) , (1) 

V(t) - eB ( ^1±^) mt) X b] = (3L^5^\ [v(0 X b] , (2) 



nil + 1712 J mi + 7712 V"^i "^2 

where qiq2j^^F {r{t)) (F = —dU/dr) is the force exerted by the particle 2 on the particle 1. The coupled, nonlinear 
differential equations (H]) and ^ completely describe the motion of the particles. For solving the scattering problem, 
they have to be integrated numerically for a complete set of the initial conditions. 

From Eqs. H]) and ([2]) follows the conservation of the parallel component of the cm velocity V(t) • b — Vo|| and the 
total energy 

(mi+m2)V'^(t) uv'^it) ,9 . . 

E = Sem + = — ^ — + + qiq2fU{v) - COUSt, (3) 

but since, in general, the relative and center of mass motions are coupled the relative i?r and cm i?cm energies are not 
conserved separately. 

In the case of two identical particles, taking here electrons, i.e. mi = m2 = m, qi = q2 = —1, the equations of 
motion considerably simplify to 

v(t)+c., [v(t)xb] - ^F(r(t)), (4) 

m 

Y{t) + Lo, [^/{t) X b] = , (5) 

with the cyclotron frequency of the electron lJc = eB /m. Here the cm- motion, Eq. ([5]), can be solved which leads to 

V [t) = Vo (t) = Foil b + Vbi [u, cos [ij^t) + [b X uj sin {uj^t)] , (6) 

where Vq±_ is the cm velocity transverse to the magnetic field direction b, and u^, — (cos (/3c, sin (/3c) is the unit 
vector perpendicular to b ((/3c is the phase of the cm transversal motion which is fixed by initial conditions). The 
velocities V^n and Vo± arc related to the particles unperturbed parallel Wqi^H ^^id transverse wqj^J- velocities [v = 1,2) 

W01||+W02|| Woi^Ui + 1;o2_lU2 

^o|| = ^ , K)±Uc = (7) 

with the unit vectors xin = (cosipu^smipu^ (ip^ are the initial phases of the particles) which fix the initial transverse 
velocities (or coordinates) of the particles 1 and 2. 

With the help of the equation of motion ([4]) and relation ^ it can be easily proven that the relative (£"1) and cm 
{Ecm) energies 

i?cm = mVlit) = m(Frij| + K,\), E, = — + fU{v) (8) 

are here conserved separately. They can be expressed by the particles initial (unperturbed) velocities Voi^ according 
to E'en! = X ("^01 ^02) , E,: = ^ (vqi — V02) . Notc that Ecm and E^- are here functions only of the difference of the 
initial phases of the particles tp = ipi — ip2. 

In the second case of BC between an electron and an heavy ion of mass M and charge Ze, i.e. qi = —1,92 = Z, 
nil = ni, m,2 = M , and assuming m/M ^ 0, /i — > m, the equations of motion ([l]) and ([2|) turn into 

v(i) + tJc [v(t) X b] = -a)c[v, X b] - ^F(r(i)), (9) 

m 

V(0 = 0, (10) 

where is the given heavy ion velocity and V(^) = = const. While the relative energy Ej- itself is not conserved 
in this case, there exists nevertheless a constant of motion 

K = — ZfU{Y) + muJcT ■ [vj y.h]^ E^ + muJcT ■ [vj x b] (11) 
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which can be easily checked with the help of Eq. ([9]). In contrast to the unmagnetized case, the relative energy transfer 
during an ion-electron collision is thus proportional to 6rj_Vi±, where Sr± and Vi± are the perpendicular components 
of the change of relative position and the ion velocity. Only for ions which move along the magnetic field direction, 
i.e. (vi_L = 0), where the magnetic term in Eq. (jlip vanishes, the relative energy is conserved. This case is very similar 
to the case of electron-electron collisions except that the cm velocity is constant. 



B. Energy loss and velocity transfer 

In the general case the rate dE^/dt at which the energy E^, = of particle v changes during the collision 

with the other particle can be obtained by multiplying the equation of motion for particle v by its velocity v,y (t) — 
V(t) + g„{iJ,/m,y)v{t), where gi = 1, g2 = — 1- As qiq2^^F (r) is the force exerted by the particle 2 on the particle 1 
the integration of this rate over the whole collision yields the energy transfer (see, [ij] for further details) 

/oo 
V(r) • F (r(r)) dr. (12) 
-oo 

assuming that for t ±oo, r{t) oo and U{r{t)) 0. According to the conservation of total energy we have 
AEi — —AE2, as it can be directly seen from Eq. p^ . 

Alternatively the energy changes can be expressed by the velocity transferred to particle 1 (or 2) during the 
collision. For this purpose we split the velocity of the vth particle, Vi,{t), into v,y(t) — Voi/(<) + Sv^{t). Here vo,y(i) 
describes the free helical motion of the particles in the magnetic field, vou ~ {Qii&B/m^) [vqu x b] = 0, while 5\'i,{t) 
is the related velocity change (with (5v^(t — )■ —00) 0). From the energy change with respect to the free motion 
5E^{t) = {m^/2){\vQ^{t) + 5v^{t)\ — VQj^(t)} the total energy transfer Ai?^ can then be expressed through the total 
velocity transfer Av^ = 5w^{t 00) = [^v{t) - ^Ovit)]t^oG and Vq^ ■ Av^ = Wo^{t) ■ 5-Vv(t)\t^oo 

AEi = ^AE2 = mi (^voi • Avi + ^Av?^ . (13) 

Employing this way of calculating the energy transfer, as e.g. in [l^, the potential U{r) has to be specified already 
at an early stage of the calculation. In Sec. IIIII we will show that Eq. allows for a more general formulation in 
which the cut-off at large distances and the regularization at small distances can be treated much easier. 

So far we considered the energy transfers of the particles in the laboratory frame. In addition the energy transfers 
AE'i and AE2 can be expressed by the change of the relative energy AE,- = A(/iv^/2) and the relative and cm 
velocity transfer A[V • v], respectively. Since AE^ = —AEcm (due to the conservation of total energy), there follows 
a relation between the energy transfers AEi and AE,- given by 

AEi = !!!l—^AE, + ^A [V • v] . (14) 
m2 + mi 

Here, the second term is the total change of the scalar product V(t) • v(i). For instance, this quantity is a constant 
for the free, unperturbed motion of identical particles where Vo(t) ■ vo(t) — (wgi — Uq2)/2 = const and is given by the 
initial velocities voi,vo2. 

In the case of identical particles the relative energy and the cm energy are conserved, i.e. AEr — —AEcm — 0, and 
thus AE'i = —AE2 = fiA [V • v]. We note however that the longitudinal {AEr\\) and transverse {AEr± = —AEj.\\) 
parts of the energy transfer AEr do not vanish, where 

AEr\\ = jAw|| (2w^|| + At-ii) ^ 0. (15) 

Here Vr\\ — foi|| — i'o2|| and Av\\ are the components of vo(t) and Av parallel to b, respectively. 

For the collision of an electron and a heavy ion, with m2 = M 3> m = mi and V = v^, the energy transfer to the 
ion AE^ = AE2 = -AEi then follows from Eq. (HH) as 

AE, = -AEr - mvi ■ Av = - (v^) - mv, ■ Av. (16) 

When the ion moves parallel to the magnetic field direction the relative energy transfer in Eq. ([TC]) vanishes and 
AE^ — — mWi||At;||, where is the component of the ion velocity parallel to b. 



5 



III. PERTURBATIVE TREATMENT. GENERAL THEORY 



A. Trajectory correction 



In this section we consider the theoretical treatment of the scattering of two identical particle, here electrons, where 
we seek an approximate solution of Eq. (jj) by assuming the interaction force between the particles as a perturbation 
to the free helical motion. For the case of ion-electron scattering the corresponding considerations and derivations 
are discussed in detail in Ref. [l^ and we thus focus on the electron-electron case in the forthcoming discussion. 

As the velocity of the cm motion is already fixed by Eq. ([6]), we have to look for the solution of Eq. (jU for the 
variables r and v in a perturbative manner 

r{t) - roW + r(i)(t) + r^^Ht)..., v(t) = Vo(i) + v^'Ht) + v^^\t)..., (17) 

where ro(<),Vo(t) are the unperturbed two-particles relative coordinate and velocity, respectively, r'^")(t), v'^"^(t) oc 
q^"F„_i {n = 1,2, ...) are the nth order perturbations of r{t) and v(t), which are proportional to g^" (for electrons 
q = — 1). Fn{t) is the nth order correction to the two-particle interaction force. Using the expansion ([T7]) for the nth 
order corrections F„ we obtain 

F (r(t)) = Fo (ro(t)) + Fi (ro(t), r,{t)) + (18) 

where 

Fo {ro{t)) = F (ro(t)) = ~i [ dk(7(k)ke'^''--»(*), (19) 



Fi(ro(i),ri(i))= [ri{t) ^ 



^) F(r) = [ dk[/(k)k [k • ri(i)] e*k ro(t)^ (20) 



r=ro(t) 

In Eqs. (|19p and (PO)) . we have introduced the two-particle interaction potential U{r) through F(r) = — (9[/(r)/(9r and 
the force corrections have been written using a Fourier transformation in space. 

We start with the zero-order unperturbed helical motion of two electrons in the relative frame 

ro{t) = Ro + hvr\\t + a [u^ sin (uJct) — [b x u^] cos {iVct)] , (21) 
Vo(t) = ro(t) , 

where = {cos tpr, sin ipr) {ipr is the initial phase of the relative motion) is the unit vector perpendicular to the 
magnetic field, Vj.\\ = t'oiH ^ ^^02[| a-nd vo±Ur — wqij^Ui — vq2±U2 (with vq± ^ 0) are the unperturbed relative velocity 
components parallel and perpendicular to b, respectively. Here Vr\\h is the relative velocity of the guiding centers of 
two particles, and a = vo±/uJc is the relative cyclotron radius. It should be noted that in Eq. (j2ip the variables u,. 
and Rq are independent and are defined by the initial conditions. Explicitly the relative cyclotron radius a and the 
phase ipr are connected to the particles cyclotron radii oi, 02 and phases ipi, ip2, and ip — ipi — ip2, according to 

n,piVl - nr,piV2 2 2 2 

a = fl]^ + flj — 2aia2 cos (22) 



a 

The equation for the first-order relative velocity correction is given by 

^^'\t)+uJ,[^r('\t) X b] = ^F(ro(t)) (23) 

m 

with the solution 

r(i)(t) ^ ^ {bP|| (t) + Re [Pi(t) - b (b • Px(i)) - z [b X P^it)]]} ■ (24) 
m 

v(i)(t) = r(i)(t). 

Here we have assumed that all corrections vanish at t —oo and have introduced the following abbreviations 

Pl\{t) ^ -i J dkU Ck)ik-h) f dT {t - t) e'^-''°^''\ (25) 

P±(t) = -— / dki7(k)k / dr[e 

t^c J J-oo 

wherein ro{t) is given by Eq. (I^T]) and the Fourier representation of the interaction force, Eq. ITOl) . has been used. 
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B. First and second order energy transfers 

The total energy change A_Ei of the first particle during collision with the particle 2 is given by Eq. (fT2)). Insertion 
of Eq. into the general expression yields 

/\Ei = A£;|^^ + A£;f ^ + (26) 

where 

/OO /* OO 

dtVo(t).F(ro(i)), AE[^^=q^f dtVo (t) ■ F, {ro{t),r,{t)) (27) 
-oo — oo 

are the first- and second order energy transfer, respectively. 

We now introduce the variable s — Ro± which is the component of Rq perpendicular to the relative velocity vector 
Vr\\h. From Eq. ((2T|) we can see that s is the distance of closest approach for the guiding centers of the two particles' 
helical motion. For practical applications the energy change is now given by the average of AEi with respect to the 
initial phases of the particles ipi, ip2 and the azimuthal angle of s. Such an averaged quantity / will be abbreviated 
by (/) in the forthcoming considerations. 

We start with calculating the first order longitudinal relative velocity transfer, Awy^^ , which also contributes to the 



second order relative energy transfer according to Eq. ()15|) . This quantity can be easily extracted from the parallel 
component of the first-order relative velocity correction, b • v'-^-'(t), in the limit t —f oo, i.e. after completion of the 
interaction. Thus, using the Fourier series of the exponential function e'^'''"^ obtain from Eqs. (p4)) and ((25)) 



Avl^'> = / dkt/(k)fc||e^'^-^" y e"(^'-^)j„(fcxa)(5(C„(k)). (28) 

II m J " f-^ 



n— — oo 



Here J„ are the Bessel functions of the nth order, CnO^) — nujc + fcjprllj = k • b and k± are the components of 
k parallel and transverse to b, respectively, tanO = ky/k^- Note that the relative cyclotron radius a as well as the 
velocity of the cm motion Vo_l depend on Lp — Lpi — (see Eq. ((T]) and ((121) )• Performing the averages with respect 
to the initial phases (pi, Lp2 and the azimuthal angle ds results in (At;^|^^) — for spherically symmetric interaction 
potentials (C/(r) = U{r) and C/(k) = U(k)) (see, e.g., Refs. [13, ll^l) due to symmetry. 

The first-order energy transfer is obtained by substituting Eqs. ([6|), (fT9|l and ([21]) into the first one of Eq. (|27p . 



Similarly as the first-order relative velocity correction A?;|'j^'' (pS)) the first-order energy change also vanishes if averaged 
with respect to the initial phases of the particles Lp2 and {}s- 

Thus the energy change receives a contribution only from higher orders, and we next evaluate the second order 
energy transfer AS^^ by inserting Eqs. (g]), ([201), (EB), (EH) and (US]) into the second equation of Eq. (gT]). This 
quantity is then averaged with respect to the initial phases of the particles ip2 and the azimuthal angle t?s of the 
impact parameter s. The obtained angular integrals are easily evaluated using the Fourier series of the exponential 

(2^ 

function. After averaging the energy transfer AE\ with respect to (pi and Lp2 the remaining part will depend on 
5 ((k + k') • b), i.e. the component of k -I- k' along the magnetic field b. Thus this (5-function enforces k -|- k' to lie in 

the plane transverse to b so that )'^° 5 ((k + k') • b) = e''^ ''5 ((k + k') • b), where Q = k_L + k'^. In addition, 

instead of phase ipi we introduce a new variable according io Lpi = ip + ip2. Then after performing ip2 -integration 
there remains average with respect to the phase difference ip. The result of the angular averaging finally reads 

(Aijf))= (29) 







where 



£{^,s) = dkdk'U{k)U{k')jQ{Qs)S{k'u+ku) y (-l)"e"(''-''') (30) 



xJn (k',a) Jn (kj_a) Gn (k, k') ( Vo||fc|| — nojc-^ cosz? ) , 

voj_ 



" ^ ' ' (c„(k') - zof (Cn(k') - *o) (C„-i(k') - zo) (c„(k') - tO) (c„+i(k') - tO) ■ ^ 
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The phase factor 'd = ipr — (fie '^i^ Eq- pO]) can be evaluated expHcitly by using Eqs. (O and (|22p . Introducing the 
phase difference of the particles if = ipi — ip2 it reads 

^tiv.-'P.) ^ ^^■9 ^ aj- al + 2iaia2 siny> ^^^^ 
-y/ (flj — a^)^ + Aafa^ sin^ (/? 

The series representation ([SO]) of the second-order energy transfer is valid for any strength of the magnetic field. 

(2) 

For most applications it is also useful to integrate the ipi,ip2, f^g-averaged energy transfer, (AE ) ) , with respect to 
the impact parameters s in the full 2D space. We thus introduce a generalized cross-section [l^, Inl, M, [13 through 
the relation 

{AE['^)sds^ Si^,s)sds. (33) 

Jo 27r 7o 2n Jq 

As CT(iy3) results from the s-integration of the energy transfer £{lp, s) (j30p one obtains an expression for ^{lp) which 
represents an infinite sum over Bessel functions. Moreover, assuming regularized interaction and performing k- 
integration in ct(</3) yields an infinite sum over modified Bessel functions (see, e.g., an example for ion-electron collision 
in Ref. ^14*1 ) . For arbitrary axially symmetric interaction potential similar expression is derived in Appendix |B] (see 
Eq. ()Bip ). However, for practical applications it is much more convenient to use an equivalent integral representation 
of the effective cross-section which does not involve any special function. This expression can be derived from the 
Bessel-function representation of the cross-section '^{f) using the integral representation of the Dirac J-function as 
well as the summation formula for '}2n^^^'^ "^n (^) H- The energy transfer a{ip) after lengthy but straightforward 
calculations then reads 

o'C'^) = o^IlM +o^-L(<y£'), (34) 



with 



^c^rll JO 



sin t\ ( t 

X ( fc^ + ^i"— ^ j sin (fc||<5t) Jo ( 2fc_Lasin - 



2 „474t/„ , r°° 



2(27r)^qV"K, 



a±{ip) = - "^""^ COS79 / tdt I dk|[/(k)|^ (36) 

TOW^|Wr||| Jo J 

sint\ , f t\ ^ f ^, , 



k']^—j— j k± COS (fc||<^i) COS ( 2 ) ( 2fc^asin 



2 

where 5 — \vr\\\/ujc is the relative pitch of the particles helices, divided by 27r. In Eq. (|34p the cross-section ct((p) has 
been splitted into two parts which correspond to the cm motion along (ct|| (ip)) and transverse {a±{ip)) to the magnetic 
field. 

An expression similar to Eqs. (|^^ and ([50]) has been obtained in [3, [23| for electron-heavy ion (no cyclotron 
motion) collision where the direction b of the magnetic field and the direction n,. = v^/vr of the relative velocity 
Vr — We||b — Vi (ugii is the component of electron velocity parallel to the magnetic field) of the electron guiding 
center is singled out in the argument of the i5-function and the summand of the n-summation. This prevents a closed 
evaluation of the energy transfer for ion-electron collision for arbitrary direction of the ion motion with respect to 
the magnetic field. But if the ion moves along the magnetic field and the guiding center Vr\\h has no component 
in transverse direction, as it is always the case for a collision of two identical gyrating particles, the energy transfer 
in such ion-electron collisions can be evaluated in the same manner as already discussed in context with deriving 
Eqs. p4p - (|36p in a straightforward manner. 

For that case of ion-electron collisions with constant v.^ = v^h the energy loss of the ion is given by AEi = 
— mti,j|| Ai;||, see Eq. p^ . because the relative energy transfer AE^ here vanishes according to Eq. pT|) . The velocity 

transfer (Av'^p) required for {AE'y^^) can be extracted from Eqs. and ([501) . Fo^' identical particles this is the term 
in Eq. ([30]) which is proportional to Vo|| (see, e.g., Eq. HH)), i.e. 

(2), _ f^"" dip 



i^')^ 2J"i(^,s) (37) 
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with 

Uiiv^s) = - ^""'f^^. I dkdk'f/(k)[/(k')A:||Jo(Qs)<5(fc[| +fc|i) (38) 

oo 

X (-l)"e™(«"«')j„(fc^a)J„(fc^a)G„(k,k'). 

n— — oo 

For the present case of these specific ion-electron coUisions the quantities and /i — m/2 in Eq. ([38|) have to be 
replaced with and m, respectively. Thus the s-integrated {Av'"^'') for ion-electron collisions corresponds to Eq. ([35|). 
i.e. more precisely to 

It should be emphasized that for ion-electron collision the phase dependent transversal relative velocity uo±(<^) is 
replaced by the electron transversal velocity Ve_L and the relative cyclotron radius vo±/uJc for identical particles is 
replaced by v^^j^c- Thus the integrands in Eqs. ()29|) and ()39|) (and similarly in other angular-averaged quantities) 
do not depend on the phase <f {£{^p, s) — S{s)) and, therefore, we have — (Ai?J ') = —S{s) — > (AiJ^ ) for the second 

(2) 

order energy transfer to the ion {AE^ ). 



C. Second order relative energy transfers 

We now turn back to the case of identical particles and consider the evaluation of the angular-averaged relative 
energy transfer in longitudinal direction, (AEf^^) = —{AE^'^^). According to Eq. ([T5|) the longitudinal relative energy 

transfer (AS^y'') involves both (Au^^'') (l37|) and the the angular-averaged value of (Aw^^'')^. The latter is calculated 
from Eq. ((28)) . The angular-averaging procedure is the same as for deriving Eqs. ((29|) and pO)) and yields 

{{Av\^^r)^ l'^^U,{^,s) (40) 

with 

U^iip.s) = "^^^l] '^^f /dkdk'C/(k)C/(k')fc^Jo(Qs)(5(fc;| +fc||) (41) 
m^ll \ J II II 

oo 

X (-l)"e"^(«-«')j„(fc^a) J„(A;»<5(C„(k)). 

n— — OO 

The relative energy change in the longitudinal direction is then calculated using Eqs. (fT5j) . ([?7|l and ([¥0)1 

(A£;f2> ^ -(Ai^^J)) =yj|'''g£.^(^,.) (42) 

with 

^ri (<^, s) = (.^, s) = -— [U2{(p, S) + 2Vr\\Ui{(p, s)] . (43) 

For further applications (see Sec. |V| also the angular averaged square of the first order relative energy transfer in the 
transversal direction, {{AE^^^)'^) is needed, which is given by 

Similarly to the definition and derivation of Eqs. (I33p and (j34p we define the related cross-sections by 

CTri_ = / {AE^;'l)sds^ / /a,^(^), (45) 

nOC /'27r J 

al, = ((A£;«)2).d. = y g^.i(^) (46) 



9 



with 



£r± {ip, s) sds 



m \ Vr\\ CO, 



2 4:14 poo 



sint 



(5t ( fc^ + fc^ J sin {k\\6t) - /c[| cos {k\\St) 



(47) 



f^rl (V) = "T^rW / U2{'f,s)sds 
4 Jo 



2 {2TTf q^^^6 dt J dk \U {k)\^ k^^ cos {ki\6t) Jo ^2k±a sin ^ 



(48) 



IV. THE ENERGY TRANSFER FOR THE SCREENED AND REGULARIZED POTENTIAL 



For the Coulomb interaction U{k) — Uc{k), the full 2D integration over the s-space results in a logarithmic 
divergence of the k-integration in Eqs. ([55)) . ([T7|) and To cure this cutoff parameters fcmin and /cmax must 

be introduced, see [13, H^l for details. But the averaged energy transfer, Eq. (|29p with Eq. ((30)) . can be evaluated 
without further approximation for any axially symmetric interaction potential, ?7(k) = J7(|fc|| |, A:j_). In this case the 
averaged energy transfer can be represented as the sum of all cyclotron harmonics as it has been done for ion-electron 
interaction in Ref . [IJ] . 

To continue the interaction must be specified. In the following we consider throughout the regularized screened 
potential U{r) — U^ir) introduced in Sec. Ill Al with 



:,(r)=(l-e-/^) 



-r/\ 



C^R(^||,fcj 



{2nf\kl+n^ kl + xV ' 



(49) 



where ^ ^2 _|_ ^-2^ ^2 ^ ^2 ^ d^^a, d^^ = X^^ + X , for which the k-integrations involved in Eqs. ((251), (gT]) 
and (|48)) converge. 



A. Second order energy transfer 

A simple but important particular case is that of vanishing cyclotron radius (woi± = vo2± = 0), i.e. when initially 
the electrons move along the magnetic field. From Eqs. (pSjl . and (PT|) . and using the regularized potential 
it is straightforward to show that in this case — and U2{f,s) ~ 0. Substituting Eq. (|49p into Eqs. ([50]) . ([55]) 

and (|43p and setting a = we obtain after k, k'-integrations 

£ [if, s) = -^£rL [v. s) = -^^^4^ [(Kis)ifi(Kis) - (xis)Xi(xis)]' , (50) 

Vr\\ mv^^^s^ 

where = S^'^ + X^^, xi = + d^^ and Kn are the modified Bessel functions. From Eq. ([50]) it is seen that the 
modified Bessel functions guarantee the convergence of the impact parameter integrated energy transfer at small s. 

Eq. (|50p is the leading term of the expansion of the the energy transfers with respect to small cyclotron radius a and 
was obtained in the limit a — *■ where the particles move initially along their guiding center trajectories. Moreover, 
for (jjc — *■ cx) also the relative pitch i5 — *■ and these trajectories are rectilinear along the lines of the magnetic field. 
The particles just pass each other along a straight line and for symmetry reasons the velocity and energy transfers 
vanish. For a finite ujc corresponding to a finite pitch the contribution of the leading term in Eq. (|50p describes the 
perturbation of the guiding centers trajectories. By the reason of symmetry the next order term must be quadratic 
~ and accounts for the finite relative cyclotron motion of the particles. 

As we discussed in Sec. IIIIBI the energy transfer ((29)) must be integrated with respect to the impact parameters s 
for practical applications. For general interaction potential this has been done in Sec. lIIIBl Eqs. (|33p - ([55|) . In general 
for a study of the convergence of the s-integrated energy transfers we note that the case with s = a is most critical 
for the convergence of the cross-sections. This is intuitively clear as the gyrating particles at |ai — 02] < s < ai + 
may hit each other on such a trajectory. This should not matter for the potential (|49)) . which has been regularized 
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near the origin for exactly that purpose. On the other hand, the energy transfer for the un-regularized potentials Uc 
and Ud diverges for s = a (see Ref. [ll] for some explicit examples). 

For the present case of the regularized and screened interaction potential Eq. ((49|l . i.e. substituting this potential 
into Eqs. ([55)1 and the impact parameter integrated cross sections are 



2qYVol\Vr 



t^dt 

WW) 



>f2 - 1 



Ti {kR (t) , K^t, t) 



(51) 



02± "01± 



t sin tdt 





-Bit) 



T3iRit),^t,t) 



-e 



T3iKR{t),KCt,t) 



4^' 

X2 - 1 

?2 _ e2j-2 , /f /„2 / \2\ „;„2 



Here ^ = 5/A = w^,, /cj^A, R^ (t) + A{a^/X^)sm\t/2), n^X/d^l + X/X, and 



(52) 



Ti {R, C, t) = 2 + 2i? - i?2 + 1 



sint 



R^ + R+1- 



i?2 



+ 3i? + 3) 



(53) 



.F2(i?,C,<) = i?+l-;^(l-^ 



/-2 

AR^ + 9R + 9- (R^ 



6i?^ + 15i?+ 15 



(54) 



(i?, C, <) = 2 + 2i? - i?2 + 1 _ 



sint 



3i? + 3) 



(55) 



T4iR,C,t) = R + l 



i?2 



sint 
~1~ 



R^ + 2R'^ + 3R + 3- (R^ + 6R^ + 15R + 15) 



(56) 



Note that the transversal cross-section a±{(p) vanish at ai = 02 and the energy transfer occurs only due to the cm 
motion along the magnetic field a\\{ip). We also note in particular the dependence of the sign of the transversal 
cross-section Eq. ([52|) on ai — a2- 

As already discussed in Sec. IIII Bl the energy transfer in the case of electron-heavy ion collisions is given by the 
s-integrated (Awl,^'), see Eq. ([39]), now with CT|[((y5) from Eq. ([51]) . 



B. Second order relative energy transfers 

Similarly, substituting Eq. ([49)) into Eqs. ([47)) and ([48)) . we obtain for the relative cross-sections CTr_L ((/?) and ct^i (</?) 



ni \vr\\ WcX Jo R'^ (t) 



dt 



yi{R{t),^t,t) 



X2 - 1 



yii^Rit),x^t,t) 



4>^2 
>f2 - 1 



J^2(>^i?(t),><^i,t) 



(57) 



dt 

W) 



-Bit) 



- 1 



V2iR(t),^t) 



+e 



Vi {xR [t) , K^t) -f — -V2 {>iR [t] , >iit) 
x'^ — 1 



(58) 



where 



V,{RX) = l-^{R + l), V2{RX)^^ 



R + i- i-: (R^ + m + i) 

R'^ 



(59) 
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yi{RX,t) = -R' + C {3R + 3 -R^)+C 1^1- 

y2{RX,t) = {C^ + l){R + l)-^{R^ + 3R + 3 

R^ + 4i?2 + 9R + 



sint 



R^ + R 



3R + 3) 



(60) 
(61) 



sint 



6R^ 



15i?+ 15 



It must be emphasized that in contrast to the quantity '^{(p) the relative cross-sections CTr±(v') and Wriiy^) do not 
depend on the cm velocity components Vo|| and Vo±. 

Next we also consider the cross-sections a{(p), ar±{(p) and Wri{(p) for vanishing cyclotron radius, a — ^ 0. In this 
limit a±{(p) = and arii^p) = as can be easily seen from Eq. as well as by explicitly evaluating Eq. with 
R{t) = ^t. The integration of Eq. (|50)) with respect to the impact parameter (see Appendix [X\ for details) yields 



a{(p) 



(Tr± W) 

tVll 



1 



(><2 _ 1) 



1 



(62) 



Note that ^((p) at a ^ can be alternatively evaluated from Eq. ([5T|) . In this limit the expression for R{t) is simplified 
to R{t) — ^t. After performing the ^-integration we again arrive at Eq. (p^ . In addition Eq. is approximately 
valid also for finite cyclotron radius a, assuming that the longitudinal velocity Vr\\ is larger than the transversal ones, 
Vo± and Vq±. Indeed, in this case R{t) ~ since 6 ^ a and in Eqs. ([5T|) and ([5^ the transversal cross-section a±{ip) 
can be neglected compared to the longitudinal one. This indicates that in the high velocity limit with v^w ^ Vo±, Vo± 
the transversal motion of the particles as well as its cm transversal motion are not important and can be neglected. 
This high velocity limit of the cross-section Wri {f ) is obtained from Eq. ((58)) . Since only the contribution of small t 
is important the function R (t) is approximated by R (t) ~ ^(^^ -|- c? ly?)^!"^ ~ ^f(l -f a? /iS"^) (here we keep the small 
term ~ a^/J^ because ari (f) vanishes at a ^ 0). Using this result for R (t) from Eqs. ([52]) . (l58|) and (f62|) in the high 
velocity limit we obtain within the leading term approximation for the cross-sections 



crri_ [(p) ^ 3 — A(>f), 



2q' 



mv, 



r\\ 



' I CU II 

r\\ 

A(>.), 



(63) 
(64) 
(65) 



^ II 
r\\ 



"HI 



It should be emphasized that the 



Here A{x) is a generalized Coulomb logarithm and is given by Eq. (|A5[) . Note that W{(p) and Wr±{(p) are isotropic, i.e. 
do not depend on ip while ari{(p) contains a term which is proportional to cos (p. Also, in the high velocity limit the 
cross-sections do not depend on the magnetic field strength and decays as CT|| (ip) ' 

and Wriiip) ~ v~^^. 

Finally we briefly turn to the case of small relative velocity, v^h ^ vo± , Vo_l ■ 
integral representations of the cross-sections, Eqs. ([5T|) . ([57|) and ([58)) . are not adopted for evaluation of these quantities 
at small velocities. For this purpose it is much more convenient to use an alternative Bessel-function representation of 
the cross-sections as shown in Appendix [Bl In addition, it is expected that the limit of small v^w is the most critical 
regime for a violation of the perturbation theory employed here. Therefore explicit analytical expressions in this 
limit can be useful for an improvement of the perturbation theory by comparing the analytical results with numerical 
simulations, see Sec. IVll 



V. IMPROVED TREATMENT FOR THE REPULSIVE INTERACTION 

The magnetic field drastically changes the scattering problem of two charged particles as discussed in Sec. ID One 
important consequence of a strong magnetic field is the backscattering of the particles from the repulsive potential 
barrier even at a finite impact parameter s (let us recall that in the case of Rutherford scattering this occurs at 
either vanishing impact parameters, s — > 0, or at vanishing relative velocity). For instance, in a strong magnetic 
field with a <C A, s and in the case of attractive interaction the velocity and energy transfers are very small and for 
symmetry reasons vanish with increasing B, see the similar discussion in Sec. IIVI However, in the case of repulsive 
interaction the magnetic field together with the interaction potential forms a potential barrier because of the particles 
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motion is effectively one-dimensional. Here two possible scattering regimes must be clearly distinguished. To this 
end, we consider an exactly solvable model for two interacting particles moving in the presence of an infinitely strong 
magnetic field on rectilinear trajectories along the field with vanishing cyclotron radii. Introducing the relative 
coordinate ( — zi{t) — Z2{t) the relative energy conservation (note that the cm energy is also conserved since V = 
bV||(i:) = bVo|| = const) for the regularized Yukawa potential can be written in the form: 



1 



so 



(l-e-'») 



-r/X 



(66) 



where <ji = \qi\/qi, ^2 = 1921/92, r'^{t) ^C,'^{t) + s'^ 



So 



2\q^q2\f/^ivl^y 



Here /i is the reduced mass, s is the impact 



parameter and Vr\\ is the initial longitudinal relative velocity of the particles. Since the relative energy is conserved 
there is no relative energy transfer and Ai?r — 0. Then the energy transfer of the particle 1 is related to the relative 
velocity transfer Au|| = C(+oo) — C(— 00) by Ai?i — /iVo||Aw||, see Eq. (fT4|). For reasons of symmetry, no velocity can 

be transferred from particle 2 to the particle 1 if the interaction is attractive (9192 < 0), see, e.g., Eq. ([5^ . This may 

1/2 

also be true for repulsive case (91(72 > 0). For instance, at so < A (or equivalently x — Vj.\\/vs > (2/1^) with v = X/X 

and — \qiq2\ ^'^ / fJ,X) and arbitrary s as well as at sq > X (or equivalently x < (2/;^)^^^) and s > Sm — Xrjmix), 
where rjm (x) is the root of the transcendental equation 



x 

Y 



(67) 



and the energy of two particles relative motion is larger than the energy of the potential barrier which again yields 
vanishing velocity and energy transfers for repulsive interaction. Thus for qiq2 > the energy transfer occurs at 
So > X and s < Sm- In this case the velocity transfer is Au|| = — 2i',.|| which corresponds to a reversion of the initial 
motion, i.e. to a backscattering event. Then the energy transfer is 



AEi = -AE2 = -2^yo||tV||0 (9192) eiv^c 



)0 (s„i - s) , 



(68) 



where Q{z) is the Heavyside function and — 2 I9192I ^"^/fJ-X. 

Consider now the energy- velocity transfers integrated with respect to the impact parameters s. The result reads 



1 

A2 



A^i(s) 



sds 



1 

A2 



AvJs) 



sds — xrf^ {x) 



(69) 



1 /2 

Here rjm is a function of x = Vr\\/vs. (Note that this function vanishes at x > (2/;^) , where the transcendental 



equation ([67]) has no solution). Consider two limiting cases. At Vj.\\/va < (2/z/)^^^ we obtain from Eq. ([67] 



Vin {x) 



3/2 



2i^ + 1 



1/2 



and the s-integrated velocity transfer vanishes as ~ — xY . At Vr\\/vs 0, Eq. ([67|) yields 



rim {x) ~ In 



2/^2 



In (2/x2) 



and the s-integrated velocity transfer vanishes like 

1 |-°°Aw||(s) 



•,ds • 



:ln^ 



2/a;2 



ln(2/a;2) 



0. 



(70) 



(71) 



(72) 



From this simple example it is clear that the perturbative treatment developed in the previous sections is not appli- 
cable for repulsive interaction and in the presence of strong magnetic field when the hard collisions like backscattering 
events may occur. Nevertheless the second order treatment for the repulsive case can be improved when the energy 
transfers due to the hard collisions are involved in the theory as the leading terms. The simple example considered 
above suggests that this can be done by exphcitly using the energy conservation, see, e.g. Eq. ([66]) . For simplicity 
consider the case when the relative energy of particles is conserved, AEr — 0. Then from Eq. (fT5|) (in the general case 
to/2 must be replaced here by the reduced mass /i) and the obvious relation A£'j,|| = —AE^.^ we obtain a quadratic 
equation for the relative velocity transfer Aw|| which has two real solutions 



Av\ 



Av = 



(73) 
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with (/) = 2Ai?rj.//it^,^j[- The velocity transfers At;J and Av^^ correspond to two different scattering regimes. In 
particular, in the limit of strong magnetic field the transversal energy transfer is small, <C 1, because of the 
transversal motion is strongly hindered and 

Avl ~ $, Avf ~ -«,|| (2 - $) , (74) 

where $ = ^(j) + ^(j)^ ■ In the limit of vanishing (j), Au^ while Av^^ — > —2vr\\ which correspond to over the barrier 
and backscattering events, respectively. Thus, in contrast to the Rutherford classical scattering the magnetic field 
may form two scattering channels, I and II. For instance, the scattering with an attractive interaction is realized in the 
regime I while in the repulsive case the scattering may occur in both regimes I and II depending on the relation between 
the initial relative energy, the height of the potential barrier and the strength of the magnetic field. In the case of the 
infinitely strong magnetic field the boundary between I and II is fixed by the arguments of the Heavyside functions 
in Eq. (p5)) . According to Eq. (1751) the velocity transfers Au|| and AvJ^ obey an exact relation Au|| + Av^^ = —2vr\\- 
We further consider the kinematically allowed range for the regimes I and II. The dimensionless transversal energy 
transfer ip can range — 0max ^0^1, where 0max — '^o±/'^r\\' ^OJ- initial transversal relative velocity. Here 

(j) = — 0max {(k — 1) corresponds to a complete transfer of the initial relative transversal (longitudinal) motion to the 
parallel (transversal) one. Therefore the quantities Au|| and Av^^ are restricted in the domains 

^ A I ^ '/'max /rrrN 

- u^ll ^ Aw|| < v^i\—===^—, (75) 

V ^ + Yraax + -L 



Vr\\ (l + V^l" 



(76) 



where the boundary between I and II is fixed by Aw|| — — 

The discussion above indicates that the second order perturbative treatment for the repulsive interaction in the 
regime II can be improved if instead of the standard approximation developed in the previous sections, the second 
relation in Eq. (|74p is used. In particular, in the presence of a strong magnetic field the modified second order velocity 
transfer reads 



A' 



,(2) 



-2Vr 



AE. 



(2) 



1 (AEi'l 



2 y fiv\ 




(77) 



Here AE"^^ and AE^^2 the first and second order relative transversal energy transfers, respectively (see, e.g., 
Eqs. (|42)) - (|44)l for collisions of two identical particles). Note that the second term in Eq. ([77)1 within the square brackets 
enters here by the opposite sign compared to the standard perturbative treatment where the first (backscattering) 
term in Eq. ([771) i^^t involved. Since we assumed the conservation of the relative energy of the particles these results 
are valid both for BC of two identical particles and for ion-electron collisions (with Vi± = 0). 



VI. COMPARISON WITH SIMULATIONS 



A. Classical trajectory Monte— Carlo (CTMC) simulations 



A fully numerical treatment is required for applications beyond the perturbative regime and for checking the validity 
of the perturbative approach outlined above. In the present case of binary ion-electron or electron-electron collisions 
in a magnetic field and with the effective interaction U^{r) (I49p the numerical evaluation of the BC energy loss 
is very complicated, but can be successfully investigated by classical trajectory Monte-Carlo (CTMC) simulations 
[Tol . [ill . [T3 |. In the CTMC method 25] the trajectories for the relative motion between the ion and an electron 
are calculated by a numerical integration of the equations of motion (i.e. Eq. ([4]) for electron-electron collisions and 
Eq. ([9]) for ion-electron collisions, respectively), starting with initial conditions for the parallel Vr\\ and the transverse 
vo± relative velocity. The initial positions are chosen to correspond to a certain impact parameter s and are located 
outside the interaction zone, which is - employing a screened interaction like (|^ - defined as a sphere of several 
screening lengths A about the ion. The numerical calculation stops after the electron has left this interaction zone, 
that is, when the collision is completed. Deducing the velocity changes from the initial and final velocities uy, v^, 
yields the energy transfer AEi (|14p (with AE,- = 0) or AE.^ ([TB)) . The required accuracy is achieved by 
using a modified Velocity- Verlet algorithm which has been specifically designed for particle propagation in a (strong) 
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r|| s r|| s 

FIG. 1: (Color online) The cross-section ai Eq. (|8ip for ion-electron (with Z > 0) collisions in terms of the dimensionless 
relative velocity component v^^^/vs at a = 0, a^/A = 0.25 (left panel) and a/A = 0.5, Os/A = 0.125 (right panel). The solid and 
dashed curves were obtained for uo — Xo/X = 5x 10~^ and I'd = 5 x 10~^, respectively. The curves with filled symbols are from 
CTMC simulations. The curves with open symbols represent the second order perturbative treatment with constant A = Aq. 
The results of the second order perturbative treatment with a modified (dynamical) cutoff parameter X{Vr\\) (see Eq. (|82|l 'l are 
given by the curves without symbols. 



magnetic field [261 127[ , and by adapting continuously the actual time-step by monitoring the constant of motion 
(HI) or K (fTTj) . The resulting relative deviations of Er or K are of the order of 10~^ — 10~^. 

The desired average over the initial phases, i.e. over the orientation of the transverse relative velocity vqj. and the 
impact parameter s is performed by a Monte-Carlo sampling [2^ [2^ of a large number of trajectories with different 
initial values. The actual number of computed trajectories is adjusted by monitoring the convergence of the averaging 
procedure. Around 10^ — 10^ trajectories are typically needed for the energy transfer for one set of initial relative and 
cm velocities and at a given magnetic field. 



B. Results 



For the forthcoming discussion we put the equation of the relative motion of two particles in a more appropriate 
dimensionless form by scaling lengths in units of the screening length A and velocities in units of a characteristic 
velocity Vg defined by 

= (78) 

Let us recall that qi ^ q2 ^ —1, fi ~ m/2 and qi ~ —1, q2 — Z, fi ^ m for electron-electron and ion-electron 
collisions, respectively. This velocity Vg gives a measure for the strength of the Coulomb interaction with respect 
to the (initial) kinetic energy of relative motion /xw^/2. For Vr < Vs the kinetic energy is small compared to the 
characteristic potential energy |gig2|^^/A in a screened Coulomb potential and we expect to be in a non-perturbative 
regime. A perturbative treatment on the other hand should be applicable for ^ v^. 

The scaled version of equations (j4]) or ([9]) only depends on the two dimensionless parameters as/A and A/A, and the 
initial conditions with positions scaled in A and velocities in Vg. Here — Vs/lOc is the cyclotron radius for v± = Vs 
and the parameter Os/A oc Vg/B represents a measure for the strength of the magnetic field compared to the strength 
of the Coulomb interaction (which is cx u^, see ([75)) ). The ratio A/ A describes the amount of softening of the screened 
interaction at r ^ with qiq2^^U-s.{r 0) ^ qxqi^ jX. 

In the analytical perturbative approach we thus apply the same scaling of length and velocities and introduce for 
electron-electron collisions the dimensionless cross-sections 
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FIG. 2; (Color online) The cross-section en Eq. (|8ip for ion-electron (with Z > 0) collision in terms of dimensionless relative 
velocity component v^^^/va at o = (left panel) and a/A = 0.1 (right panel). The curves with and without symbols correspond to 
CTMC simulations and the second order perturbative treatment, respectively. Left panel, Os/A = 0.5 (solid curve), Os/A — 0.25 
(dashed curve), a^/A — 0.125 (dotted curve). Right panel, fls/A — 0.1 (solid curve), a^/A = 0.05 (dashed curve), Os/A = 0.025 
(dotted curve). 




V ,,/v 

r|| s 



FIG. 3: (Color online) Same as in Fig. [2]but for a/A — 0.25 (left panel) and a/A = 0.5 (right panel). In left and right panels, 
as/\ = 0.5 (solid), a^/A = 0.25 (dashed), a^/A = 0.125 (dotted). 



The corresponding dimensionless cross-section for ion-electron interaction is given by 

a, = / {AE,)sds^ -2/ (Aw||)sds. (81) 

TOWij|WsA^ Jo WsA^ Jo 

It should be noted that in the regimes where the cross-section CT|| (ip) is not strongly sensitive with respect to the 
initial phases ip the quantity ai in Eq. (jSip is approximately given by ai ~ 2oji(see also the relation (1391) ) in the units 
introduced above. An example of such regime has been considered in Sec. lIVBl where CT|| {(p) in the high- velocity limit 
and in leading order is independent on (p, see Eq. ([M]). One can expect that the approximate relation between the 
ion-electron and electron-electron cross-sections is valid within second order perturbative treatment with arbitrary Z 
and within exact CTMC simulations, where, however, the ion-electron interaction must be repulsive {Z < 0). 

Next we specify the cutoff parameter A which is a measure of softening of the interaction potential at short distances. 
As we discussed in previous sections the regularization in the potential is sufhcient to guarantee the existence of 
the s-integrated energy transfers, see e.g. Eq. ([33)l . but there remains the problem of treating hard collisions. For a 
perturbation treatment the change in relative velocity must be small compared to w,. and this condition is increasingly 
difficult to fulfill in the regime Vr 0. This suggests a physically reasonable procedure: The potential must be 
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FIG. 4: (Color online) The cross-section ai Eq. (|8ip for ion-electron collision with Z < a,t a — (left panel) and a/A = 0.1 
(right panel). The curves with symbols correspond to CTMC simulations. Left panel, as/X = 0.5 (filled circles), Os/A = 0.25 
(squares), as/A — 0.125 (open circles). Right panel, as/A = 0.1 (filled circles), as/A — 0.05 (squares), as/A = 0.025 (open 
circles). The dashed curves are obtained employing Eq. (|69p with constant cutoff parameter v — vo ~ 5 x 10^^. 
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FIG. 5: (Color online) Same as in Fig. |4]but for a/A — 0.25 (left panel) and a/A = 0.5 (right panel). In left and right panels, 
as/A — 0.5 (filled circles), as/A = 0.25 (squares), as/X — 0.125 (open circles). 



softened near the origin. In fact the parameter A which describes the effects of quantum diffraction should be related 
to the de Broglie wavelength which is inversely proportional to Vr- Here within classical picture of collisions we employ 
in a perturbative treatment the dynamical cutoff parameter >c(vr\\) = 1 + X/X(vr\\), where 

2 2 2 l'i'l'Z2 ll^'^ 

X (t;,.ji) = C6oK-|i) + Aq, 6oK|i) = „/„2 , „2 v (82) 

In the case of ion-electron collision |(7ig2|, fJ- and vq± in Eq. (|82|) have to be replaced by \Z\, m and Ve±, respectively. 
Here Xq is some constant cutoff parameter, and 6o(''^r||) is the distance of closest approach of two charged particles in the 
absence of a magnetic field. Also in Eq. (|82p we have introduced an additional fitting parameter C. For determining 
C we consider the second order transport cross-section trtr = —{2tt / iiVr\\VQ\\)a , where a is given by Eq. (j33p . For 
the regularized potential and in high- velocity limit a is given by Eq. (j64p . where the generalized Coulomb logarithm 
A{>i) and the cutoff x = x{vr\\) are determined by Eqs. (|A5p and respectively. Setting Aq = the obtained 

high- velocity transport cross-section is then compared with an exact asymptotic expression derived in Ref. [30[ for the 
Yukawa-type (i.e. with A — > 0) interaction potential which yields C = e'^'^~^ /A ~ 0.292, where 7 is Euler's constant. 
As will be shown below the second order cross-sections with dynamical cutoff parameter (|82p excellently agree with 
CTMC simulations at high velocities. The CTMC simulations have been carried out with constant A = Aq <C A 
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FIG. 6: (Color online) Comparison of the CTMC simulations for the attractive (circles) and repulsive (squares) ion-electron 
interactions at Os/A — 0.25. The solid curves correspond to the second order theory with modified cutoff parameter. In left 
and right panels a — and a/ A — 0.5, respectively. 



a,/?i = 0.1,a, = 

aJX = 0.05 
a a = 0.025 




FIG. 7: (Color online) The cross-section Eq. (|79p for electron-electron collision in terms of dimensionless relative velocity 
component Vr\\/vs for a2 = and at ai = (left panel) and ai/A = 0.1 (right panel). The curves with and without symbols 
correspond to CTMC simulations and the second order perturbative treatment, respectively. Left panel, a^/A = 0.5 (solid), 
as/\ = 0.25 (dashed), a^/X = 0.125 (dotted). Right panel, a^/A = 0.1 (solid), a^/A = 0.05 (dashed), a^/A = 0.025 (dotted). 



(or X = >fo — -^/^o ^ 1), that is, the interaction is almost Coulomb at short distances. As an example in Fig. [T] 
we compare the cross-sections <Ji obtained with CTMC simulations (curves with filled symbols) and within second 
order perturbative treatment either with constant A = Aq (open symbols) or dynamical cutoff parameters (|82p (curves 
without symbols). Two distinct cases are considered here with vq = ^o/A = 5 x 10^^ and I'o = 5 x 10^^. Also 
a = (no cyclotron motion at the initial state), Og/A — 0.25 and a/A — 0.5, a^/A — 0.125 in left and right panels, 
respectively. It is seen that with respect to the 'softness' of the regularized potential the agreement between CTMC 
and the perturbative treatment improves, keeping all other parameters fixed, for increasing A/A, that is a weaker 
interaction potential C/r. 

Systematic investigations and comparisons of the cross-sections determined by the CTMC simulations and the 
second order perturbative treatment Eqs. ((5T|l . (|52|) . ([79|) - ((8T|) with the dynamical cutoff parameter X{vr\\) (|82|l are 
presented in Figs. IMTTI We have also compared the CTMC results with the simple (exact) model given by Eq. ((69|) 
for a repulsive ion-electron interactions, Figs. [3] and [S] Shown are ai for attractive (with Z > 0) and repulsive (with 
Z < 0) interactions (Figs.HHS]), o-\\ (Figs-jTHHl), (Figs.[TUland [TT|) as functions of Vr\\/vs for fixed cyclotron radii 

a, ai and 02 and varying the strength of the magnetic field Os/A = Bg/B with Bs — mVs/eX. For each pair of fixed Oj 
(i — 1,2) and the transversal relative velocity is determined by voi±/vs = ai/a^ (or Ve±/vs = a/as for ion-electron 
collisions). The CTMC results are indicated by the curves with filled or open symbols and the corresponding second 
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V ,,/v V ,,/v 

r|| s r|| s 

FIG. 8; (Color online) Same as in Fig. [7]but for ai/A — 0.25 (left panel) and ai/A — 0.5 (right panel). In left and right panels, 
as/\ = 0.5 (solid), a^/A = 0.25 (dashed), a^/A = 0.125 (dotted). 




FIG. 9: (Color online) Same as in Fig. [7] but for 02 — ai/2 with ai/A = 0.05 (left panel) and ai/X — 0.5 (right panel). Left 
panel, as/A = 0.05 (solid curve), a^/X — 0.025 (dashed curve), a^/X — 0.0125 (dotted curve). Right panel, a^/X — 0.5 (solid), 
Qs/X = 0.25 (dashed), a^/A = 0.125 (dotted). 

order predictions are given by the curves without symbols. Both the CTMC and second order calculations have been 
done for a regularized potential [/r with i^q — 5 x 10^^. Note also the logarithmic scales for the cross-sections in 
Figs.mHandEHni 

First we discuss some general observations of the behavior of the cross-sections. In Figs. [U [3] and [BHTT] we can 
clearly observe that in the regimes of large relative velocities Vr\\/vs > k (where typically 2 ^ k ^ 4) the second 
order perturbative treatment agrees perfectly with the numerical CTMC results. In addition in the limit of very large 
velocities Vr\\/vg 3> 1 the cross-sections (Ji and (T|| calculated either within perturbation theory or CTMC method 
with different strength of the magnetic field and transversal velocities converge to the same value. This behavior 
agrees with the predictions of the asymptotic expression l|64p which is independent on B and uoi±: vo2± (or Uej_). It 
can also be seen that the smaller the transversal velocities, the better is the convergence to the regime of Eq. (|64p . At 
large relative velocities the second order and CTMC cross-sections a± shown in Figs. [TOl and [Til agree with Eq. (|80|) 
with the asymptotic expression (|65p . Since at Vr\\/vs S> 1 the quantity a± is not affected by the magnetic field but 
behaves as a± ~ '^oi±j ^o2± ^"-"^ fixed electron cyclotron radii it will be larger for smaller (more precisely a± aj^) 
as shown in Figs. [TO] andfTTl 

At small velocities with Vr^^/vg < 1 the second order treatment considerably deviates from CTMC simulations, see 
Figs. [THSl and Figs. [SHTT] Here the second order cross-sections are given by approximate expressions (|B6|) and (|B12[) 
where the parameter x at small relative velocities is given by ><• ~ ^(0) = 1-1- A/A(0) and A(0) is the dynamical 
cutoff X{vr\\) Eq. ([82]) at Vr\\ — 0. Note that at finite cyclotron radii of the particles the quantity A(0) is a constant 
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FIG. 10: (Color online) The cross-section a± Eq. (|80p for electron-electron collision in terms of dimensionless relative velocity 
component v^\\/vs for a2 = and at ai/A = 0.25 (left panel) and ai/A = 0.5 (right panel). The curves with and without 
symbols correspond to CTMC simulations and the second order perturbative treatment, respectively. In left and right panels, 
aa/\ — 0.5 (solid curve), a^/A = 0.25 (dashed), a^/A = 0.125 (dotted). 




FIG. 11: (Color online) Same as in Fig. llOl but for 02 — ai/2 and at ai/A = 0.05 (left panel) and ai/A = 0.5 (right panel). Left 
panel, a^/A — 0.05 (solid curve), a^/A — 0.025 (dashed), a^/A — 0.0125 (dotted). Right panel, as/A = 0.5 (solid), a^/A = 0.25 
(dashed), a^/A = 0.125 (dotted). 

depending on the value of i/q and the transversal velocities. However, for vanishing cyclotron radii (as, e.g. in the left 
panels of Figs. [U [H [5] and [7]) the cutoff parameter ([5^ at small velocities behaves as A/A ~ {vs/vr\\)'^- The quantity 
{>c— 1)^ involved in Eqs. (|B6p and (|B12p falls as ^ (w^n/ws)'*. This results in a strong self-cutting at small velocities. 
Thus employing the cutoff (|82p the second order cross-section cr|| is strongly reduced and decreases as a\\ ~ f and 
(T|| ~ ii^ii at a = and a 7^ 0, respectively. The transversal cross-section <7±, i.e. the second term in Eq. (jB6|l . does 
not contain a term (x: — 1)^ and diverges as a± ~ ^r\\^ Figs. [TOl and fTTl In this small velocity regime the second 
order perturbative treatment is clearly invalid and a non-perturbative description is required. 

In Figs. [4] and [5] we demonstrate the cross-section ai for repulsive {Z < 0) ion-electron interaction. The dashed 
curves represent the cross-sections obtained from simple model considered in Sec.|Vl see Eq. (|69p . with constant cutoff 
parameter i/q. Let us recall that the model completely ignores the cyclotron motion of the particles. It is seen 
that in Figs. [4] and [5] (left panel) the agreement with CTMC simulations is quite satisfactory even for finite cyclotron 
radius a and magnetic field. However, with increasing a the CTMC simulations show a more involved picture as e.g. 
in Fig. [5] (right panel) than the predictions of a simple model (I69|) . And the deviations from CTMC becomes more 
pronounced with increasing electron transversal velocity Ve±_ and magnetic field, see e.g. open symbols in the right 
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panel of Fig. [5l In the CTMC simulations the cross-section shrinks strongly at tvji > Vs with increasing relative 
velocity. This feature can be explain on the basis of a simple model discussed in Sec.|Vl At strong but finite magnetic 
field the hard collisions like backscattering events may also occur but the total relative energy Atv^||/2 in Eq. ([66]) will 
be replaced here by the total energy +Ugj^)/2. This will reduce the domain of the backscattering events (regime 
II introduced in Sec.|V|. With increasing Vr\\ this domain will be further shrunk and finally the scattering may occur 
only in the regime I where a strong magnetic field may strongly reduce the energy transfer. (Let us recall that the ion 
moves along the magnetic field. In this case the energy transfer vanishes with increasing B for symmetry reason, see 
Refs. [13, HSl)- Obviously this effect is pronounced by increasing transversal velocity as shown in Fig.O The similar 
feature is observed also for electron-electron interactions shown in the right panels of Figs. [8] and [9l 

For ion-electron collisions we also compare the cross-sections for attractive {Z > 0) and repulsive {Z < 0) interac- 
tions in Fig.[Sl The solid curves represent the second order energy transfer which is quadratic in Z. For small velocities 
the cross-section (t.; for repulsive interaction considerably exceeds the cross-section for attractive interaction. This is 
because of the backscattering events at Z < with large energy-velocity transfers. However, with increasing relative 
velocity the agreement between second order theory and CTMC and between two CTMC with positive and negative 
Z continuously improves and is almost perfect over the full range of the velocity at v^w/vs ^ 1. 

As has been emphasized in the preceding sections there are similarities between electron-electron and ion-electron 
(with V,; — Wi||b) interactions but the energy transfers are not quite the same. In particular, the approximate relation 
ai ~ 2ct|| for the negative ions with Z < may be violated at small velocities where the electron-electron cross- 
section {(fi) is not in general isotropic with respect to (f. In Figs. [7]and[8]the cyclotron radius of the second electron 
vanishes a2 = and thus a = ai as well as o'||(</3) are isotropic. In this case the relation between the cross-sections 
is exact ai = 2a\\. In Fig. [9] we show an anisotropic case where a2 = ai/2 and ai varies from a small (left panel) 
to a large value (right panel). Another feature which is absent in ion-electron collisions is the cm cyclotron motion 
for two-electron collisions which causes the transversal energy transfer a± shown in Figs. [10] and [TT] For symmetry 
reason this quantity vanishes at ai — 02 both in CTMC simulations (within the unavoidable numerical fluctuations) 
and the second order theory, see Eq. ([5^ . In the examples shown in Figs. [TUl and [TT] the cyclotron radius 02 of the 
second electron is smaller ai , 02 < ai . Since the cross-section (t± is positive the transverse energy of the first electron 
with oi > 02 transfers to the energy of parallel and transverse motion of the second electron. 

VII. CONCLUSION 

In this paper we have investigated the binary collisions (BC) of two charged particles in the presence of constant 
magnetic field employing second order perturbation theory and classical trajectory Monte-Carlo (CTMC) simulations. 
Two distinct cases with symmetric and strongly asymmetric charges and masses of the particles have been considered 
in detail: (i) BC between two identical particles (e.g. electrons) and (ii) between an electron and heavy ion which 
moves with rectilinear trajectory (no cyclotron motion) along the magnetic field. Physically these two distinct cases 
are similar except the time-dependent center of mass cyclotron motion in the case of two identical particles. The 
second-order energy transfers for two-particles collision is calculated with the help of an improved BC treatment which 
is valid for any strength of the magnetic field and thus involves all cyclotron harmonics of the particles motion. For 
further applications (e.g., in cooling of ion beams, transport phenomena in magnetized plasmas) the actual calculations 
of the energy transfers have been done with a screened interaction potential which is regularized at the origin. The 
use of that potential can be viewed as an alternative to the standard cutoff procedure. For a repulsive ion-electron 
collisions we furthermore presented an exact solution for the energy transfer in the presence of an infinitely strong 
magnetic field. Here backscattering events may occur when the relative velocity transfer is independent on the strength 
of the Coulomb (regularized and screened) interaction Avy = —2vr\\, where Vr\\ is the initial relative velocity of the 
guiding centers of the particles. The second order perturbative treatment is clearly invalid in this case. It has been 
shown that for repulsive interaction of the particles in a strong but finite magnetic field two scattering regimes must 
be distinguished with and without backscattering events with large energy-velocity transfers. This also suggests an 
improved perturbative treatment for repulsive interaction and in the case of strong magnetic field. 

For checking the validity of the perturbative approach and also for applications beyond the perturbative regime 
we have employed numerical CTMC simulations. These CTMC calculations have been performed in a wide range 
of parameters (magnetic field and the relative velocities of the particles) and for a small regularization parameter, 
that is, for an interaction which is rather close to Coulomb at short distances. Within the second order treatment 
we have introduced a dynamical cutoff parameter which substantially improves the agreement of the theory with 
CTMC simulations. From a comparison with the non-perturbative CTMC simulations we have found as a quite 
general rule which is widely independent of the magnetic field strength that the predictions of the second order 
perturbative treatment are very accurate for v^^^/vg ^ 4 for all studied parameters and cases, with the characteristic 
velocity Vs given by Eq. (j78p . In contrast, for low relative velocities Vr\\/vs ^ 1 the results obtained from perturbation 
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theory strongly deviate from the CTMC simulations. Moreover in this regime the CTMC calculations display a large 
difference between positive and negative ions which disappears for high velocities Vr\\. Such a difference is completely 
absent in the second order perturbative treatment where the energy transfers is quadratic in the ion charge Ze. We 
have also tested the exact analytical model Eq. (j69p derived for repulsive interaction and an infinitely strong magnetic 
field by comparing it in Figs. [4] and [5] with the CTMC simulations and found that the agreement is rather satisfactory 
even for finite (but strong) magnetic fields. 

We believe that our theoretical findings will be useful for the interpretation of experimental investigations. Here, it 
is of particular interest to study some macroscopic physical quantities on the basis of the presented theoretical model 
such as cooling forces in storage rings and traps, stopping power of ion beams as well as transport coefficients in 
strongly magnetized plasmas. These studies require an average of the energy or velocity transfers with respect to the 
velocity distribution of the electrons. The cooling forces obtained by the perturbative approach are expected to be 
quite accurate if the low velocity regime only slightly contributes to the v^-average over (AE). That is, if the typical 
Vr\\, given by the maximum of the thermal electron velocity and the ion velocity, are large compared to Vs, as it is 
usually the case for e.g. electron cooling in storage rings. 

Another interesting issue not considered here is the interaction of the magnetized electrons with light ions (in 
particular with positrons, protons and antiprotons) when the magnetic field is so strong that the cyclotron motion of 
the ion cannot be neglected anymore. In this case the relative and center of mass motions are coupled to each other 
and the center of mass velocity of the particles cannot be represented in the form of simple cyclotron motion, Eq. 
A detailed comparison with cooling force measurements and a study of other aspects are in progress and the results 
will be reported elsewhere. 
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APPENDIX A: DERIVATION OF THE COULOMB LOGARITHM A{i 



In this Appendix we evaluate the effective cross section in the case of vanishing relative cyclotron radius. Using 
the definition of a{ip), Eq. (|33p . and the energy transfer Eq. ([50]) we obtain 



cr(^) = -- 



A(u) 



(Al) 



Here w = Xi//s;i and we have introduced the generahzed Coulomb logarithm 

A(u) = / [Ki (s) - uKi {us)f sds. 
Jo 



(A2) 



The integration in Eq. (|A2|) can be done using the indefinite integrals of the functions sKi{s), sK1{us) and 
sKi{s)Ki{us) [24]. This yields 
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- 1. 



Since at s ^ the modified Bessel functions behave as (see, e.g., Q) K2 (s) = 2/ s"^ ~\- 0{\), Ki (s) = l/s + 0(slns), 
Kq (s) = ln(2/s) + 0(1), the expression in the large brackets in the last line of Eq. (|A3|) vanishes as in^ g ^ 0. 
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There remains 



A (u) = [Ko (s) - Ko - 1. 

Since Kq (s) = ln(2/s) + 0(1) at s from Eq. (K^ we finally obtain 

2 + 1 



A(u) 
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APPENDIX B: THE ENERGY TRANSFER IN A SMALL VELOCITY LIMIT 



For the second order BC treatment the most critical situation is the small velocity regime where we expect some 
deviations from the non-perturbative CTMC simulations. For the improvement of the theoretical approach it is 
therefore imperative to investigate the energy transfer in the small velocity limit, \vr\\ \ ^ wo± or alternatively S ^ a. 
In principle this limit can be evaluated using the integral representation of the cross-sections, Eqs. ([SI]) . ([52]) . ((57|) 
and (|58p . However, while these expressions are very convenient to calculate the high velocity limit of the energy 
transfers (see Sec. IIV[) they are not adopted for the evaluation of the small velocity limit due to the oscillatory nature 
of the function R(t) at Vj.\\ 0. In this Appendix we consider instead an alternative but equivalent expressions 
for the effective cross-sections. For the axially symmetric interaction potential the effective cross-sections can be 
evaluated using Eqs. dSni), (ESI), dm) and dlSl). We refer the reader to Refs. [H, iO] for details. The integration of 
these expressions with respect to the impact parameter s yields the two-dimensional (5-function, <5(Q) = S(\i± + k'j^). 
Combining this function with the one-dimensional (5-function S{k\\ + fc||) in Eqs. (I30|) . ([38|) and (|4T|) yields a three- 
dimensional 5-function 5(k-|-k') and the k'-integration in the energy transfers can be performed exactly. Furthermore 
it can be shown that only the imaginary part Im[G'„(k, — k)] contributes to the cross-sections. From Eq. ([31]) follows 
that the imaginary part of G„(k, — k) is expressed by the Dirac functions 5(Cn(k)) which allows to perform the 
fc|| -integration. The final result reads 
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a.i (^) = ^ ^ [n2$„^, (fell , a)\ , (B3) 



where the function <i>n;m(fc|| , a) is defined as 



n=l 



^4 ,.00 



$„;m(fc||,a) = / V {k\\,k^) Jl{k^a)Vldk^. (B4) 



In Eq. (jBip we have introduced the notations r\n — \ — ^(5„o, = Vpn — {yr\\ /2)./(y) with /((ys) = (af — a^la}[(p). 

Equations (jBip - (|B3|) are equivalent to the integral representations ([34)1 - (|36|) . ((47l) and (|48| of the cross-sections, 
respectively. The function (fc|| , a) is taken at fc|| —n/5 which in the limit of small velocities becomes very large 
(except the term with n = in Eq. (|Bip). For the regularized interaction potential from Eq. (|B4|) at fcya S> 1 
we obtain 

<i,„.^(fc 1^^— ^, $„.3(fc|, a) - — ^""'"l^ . (B5) 
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Substituting these expressions into Eqs. (|Bip - (|B2p in the lowest order with respect to Vr\\ (or S) we arrive at 
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Here C(z) is the Riemann function with ^(5) ~ 1.0369. The function V'(w) is expressed by the modified Bessel functions 
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Note two hmiting cases of this function. At u <C 1 (vanishing cyclotron radius) and >c ^ 1 from Eq. (jB9P we obtain 
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respectively. The case of ion-electron BC is easily recovered from the obtained expressions as described at the end of 
Sec. IIVI Note that the cm cyclotron motion is absent here and in Eqs. (|Bip and (jB6p the terms proportional to the 
function f{ip) must be neglected. 

From Eqs. (jBG^ - IIBSP we then obtain that the relative transversal cross-sections at small relative velocities behave 
as Ur^_('p) ~ iij?!!, 'Wri{f) ~ f^||. The first term in the cross-section 'W{ip) vanishes as ^ while the second one which 

corresponds to the zero harmonic with n = in Eq. (|B1[) diverges as ~ This later term vanishes at Oi = 02 and 
predicts an infinitely large energy transfer at ai > a2 (or energy gain at ai < a2). As expected the cross-sections 
(|B6P - (|B8P are strongly anisotropic with respect to the phase (/s. 

Finally, we consider the small velocity limit of the cross-sections a {ip) and 'Wri_ {^p) at vanishing cyclotron radius, 
a = (see Eq. 



(7 {(p) = —(7r± (P) 



(B12) 



In this case the cross-sections behave as a{Lp) ~ Vr\\ and ar^ip) ~ w^y at small relative velocity Vr\\ (cf. with Eqs. (|B6 
and ([BT])). 
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